
import matplotlib.pyplot as plt

from GenomeReader import readFastq
from tools import naive
from tools import reverseComplement
from tools import naive_2mm
from tools import findGcByPos


seqs , quals = readFastq("data/ERR037900_1.first1000.fastq")

gc = findGcByPos(seqs[:])

for i in range(len(gc)):
    if gc[i] < 0.3:
        print(i, seqs[i])

plt.plot(range(len(gc)), gc)
plt.show()





